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The fabrication of high-quality thin superconducting films is essential for single-photon detectors. Their 
device performance is crucially affected by their material parameters, thus requiring reliable and nondestructive 
characterization methods after the fabrication and patterning processes. Important material parameters to know 
are the resistivity, superconducting transition temperature, relaxation time of quasiparticles, and uniformity of 
patterned wires. In this work, we characterize thin NbN films by using transport measurements in magnetic 
fields. We show that from the instability of vortex motion at high currents in the flux-flow state of the IV 
characteristic, the inelastic life time of quasiparticles can be determined to be about 2 ns. Additionally, from the 
depinning transition of vortices at low currents, as a function of magnetic field, the size distribution of grains 
can be extracted. This size distribution is found to be in agreement with the film morphology obtained from 
scanning electron microscopy and high-resolution transmission electron microscopy images. 

PACS numbers: 85.25.Oj, 74.25.Wx, 74.40.Gh 



I. INTRODUCTION 

The potential applications of the single-photon detector 
(SPD) in quantum cryptography, ultrafast photon detection 
experiments, and dark sky observations are very promising. 
In particular, superconductor-based SPDs have attracted con- 
siderable attention in the last decade.Wi^ Among these, re- 
search focused on the NbN superconducting nanowire SPD 
(SNSPD) for the following reasons: 1) the superconducting 
energy gap, A, is two orders of magnitude smaller than in 
semiconductor-based SPDs, which allows for the detection of 
low energy photons in the far-infrared region of the spectrum; 
2) the energy relaxation time between electrons and phonons 
is fast, approximately 17 ps, which means small delay times 
for reset with gigahertz count rates, 1 3) low dark count rates 
are attainabl e, becau se the SNSPD is operated in a cryogenic 
environment! 5 ! 6 ! 15 ! 1 ^ So far, most of the work has focused on 
the NbN SNSPD, because thin NbN has an extremely short 
superconducting coherence length of a few nanometers, £ ~ 4 
run, with a relatively high superconducting transition temper- 
ature, T c ~ 14 K, and strong electron-phonon coupling for 
fast energy relaxation times. The small £ permits one to re- 
duce the dimensionality of SNSPDs to nanoscale-sized wires 
for increased sensitivity to infrared photons with wavelengths 
A > 1500 nm. Since the superconducting condensation energy 
density per unit volume is materials specific and allows for lit- 
tle variability, given the constraints listed above, the obvious 
dimensional tunability of the device is to reduce the volume 
element of the detector that needs to go normal to trigger a 
photon count. Recently, the sensitivity of SNSPDs to ions at 
low energy and soft x-rays was explored EES p or a review of 
SNSPDs based on other superconductors see e.g. Ref. Q4] 

The operating principle of the NbN SNSPD is as 
The nanowire is biased by a DC current close 
to the critical current. When an incident photon interacts with 
the NbN nanowire, it excites a cloud of quasiparticles, that 
diffuses and drives a belt-like normal region across the wire. 
When this extended normal region appears, the resistance of 



the NbN nanowire becomes much larger than that of a paral- 
lel shunt resistor. As a result, the current redistributes to the 
shunt and a voltage pulse is detected. This process is very fast 
and takes place within ~ 10 ps. In the meantime the normal 
region in the nanowire recovers back to the superconducting 
state and the bias current flows back into the nanowire. The 
redistribution of the current at this stage is slow (1-10 ns) due 
to the large kinetic inductance of the NbN nanowire and the 
shunt resistance. During this time the nanowire cools down 
to the bath temperature (with phonon escape time ~ 160 ps) 
and the SNSPD is again ready for the detection of incoming 
photons. In the absence of incident photons, some part of 
the nanowire may become normal as well because of ther- 
mal fluctuations, which cause the detection of so-called dark 
counts ! 6 l 15 l 20 llt was proposed that the dominant contribution to 
the dark count rate is from the crossing of single vortices in the 
NbN nanowire due to thermal fluctuations in the metastable 
DC-biased superconducting state r 1 ! 22 ! 

The operation of the NbN SNSPD involves the excita- 
tion and relaxation of quasiparticles, a complex nonequilib- 
rium problem. Thus the inelastic lifetime t of quasiparticles 
plays an important role in determining the physically limiting 
SNSPD performance. On the other hand, grain boundaries are 
inevitably introduced during the growth process of thin NbN 
films. These boundaries work as pinning centers for vortices, 
thus affecting the vortex crossing in nanowires and their dark 
count rate. We expect that knowledge of the inelastic relax- 
ation time of quasiparticles and the size distribution of grains 
in thin films are important for the device optimization of NbN 
SNSPDs. 

In this work, we measure the standard materials properties 
given by the normal-state resistivity p(T) and superconduct- 
ing transition temperature T c . These are supplemented by the 
extraction of the inelastic relaxation time r of quasiparticles 
and the distribution of grain sizes P(L) in thin NbN super- 
conductors from transport measurements of the IV character- 
istics. The former is extracted from the instability of vortex 
motion in the flux-flow state, while the latter is obtained from 
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the depinning transition of vortices. The inelastic relaxation 
time of quasiparticles in the vortex state is about t ~ 2 ns in 
our NbN films. The nondestructive determination of the grain 
size distribution is dominated by domains of linear dimension 
of less than 5£ with an exponential tail for domains larger than 
~ 50£. This result is confirmed by morphology studies of the 
NbN film with electron microscopy. 

The remaining part of the paper is organized as follows. 
Sec. II describes thin film growth, fabrication, micropattern- 
ing, morphology analysis, and standard transport characteri- 
zation. In Sec. III. A, we discuss the instability of the flux- 
flow state and inelastic quasiparticle relaxation time. This is 
followed by Sec. HLB, with the investigation of the depinning 
transition and nondestructive extraction of grain size distribu- 
tion. The paper concludes with a short summary in Sec. IV. 



II. THIN FILM GROWTH AND CHARACTERIZATION 

A. Growth and fabrication 

The thin NbN films in this study were grown by a 
molecular-beam-epitaxy-type growth process called Energetic 
Neutral Atom Lithography and Epitaxy (ENABLE). EN- 
ABLE utilizes an energetic beam of neutral N atoms (kinetic 
energies of 1 to 5 eV) to activate nitride thin-film growth. 
The high energy and reactivity of N atoms allow for growth 
of high-quality, uniform crystalline thin films that are dif- 
ficult to grow by conventional chemical vapor deposition 
(CVD) or molecular beam epitaxy (MBE) techniques. EN- 
ABLE has previously been used for the growth of Ill-nitride 
semiconductors. 25 27 The thin film was grown on a 50 mm 
diameter substrate of c-axis sapphire wafer that was first pre- 
nitrided at 400° C under the ENABLE N-atom beam for 20 
minutes. Following pre-nitridation, the NbN film was grown 
at 600° C by having the Nb metal flux and N atom beams con- 
currently bombard the substrate. The Nb flux was provided by 
a Nb rod in a miniature electron-beam cell manufactured by 
Mantis Deposition Ltd (Oxfordshire, UK). The resulting NbN 
film was cooled to ambient temperature in vacuum and was 
d = 11.2 nm thick. Finally, the fabricated NbN micron-sized 
wires were photolithographically defined using MicroChemi- 
cals AZ 5214E photoresist and Microposit MF 319 developer. 
The pattern defined in the photoresist was transferred into the 
underlying film using a 50% chlorine in argon ICP etch. The 
contact pads were patterned via a lift-off technique with the 
metal (5 nm Ti, 200 nm Au) deposited by an electron-beam 
evaporator. The wafer was then diced using a resin bonded 
diamond blade. 



B. Film structure and morphology 

Structural characterization of the films was accomplished 
with scanning electron microscopy (SEM) and transmission 
and scanning transmission electron microscopy (TEM and 
STEM). The film morphology of ENABLE-grown films is 
shown in Fig. [T] The SEM image in panel (a) shows the 
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FIG. 1. The morphology of ENABLE-grown NbN film is visualized 
in (a) with the lxl yum 2 plane-view SEM image. The film thickness 
of d = 11 .2+0.4 nm and grain size and film morphology are shown in 
(b) with the bright-field TEM cross-sectional image. Bending of the 
lattice fringes around surface irregularities is shown by the Z-contrast 
STEM imaging in (c). The HTEM image of a sister sample in (d) 
shows the atomic structure of the interface and grain boundaries in 
these NbN films. 



uniformity of the film in a l//m x l//m view. The contrast 
giving rise to the grainy morphology arises from the 1 to 2 
unit-cell deep grooves above the grain boundaries in the film. 
This grain boundary morphology is shown in the TEM cross- 
sectional view of this film in Fig.[TJb). The average thickness 
of the film from the TEM measurements was d = 1 1 .2 + 0.4 
nm. The image was taken near a zone axis. Some of the 
NbN grains are crystallographically aligned to the beam to 
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FIG. 2. (color online). Grain-size analysis of the SEM image shown 
in Fig.[TJa). Panel (a): Post-processed SEM image where gray scales 
have been mapped onto a binary black and white map. Black do- 
mains are delineated by white "background" for follow-on analysis. 
Panel (b): The "Analyze Particles" method in ImageJ^^ was used 
to create the histogram of the grain-size distribution of the black do- 
mains in the lxl jim 2 area shown in (a). For small domains the 
distribution obeys a ~ 1/L 3 behavior (red solid line) with minimal 
linear dimension of L,„,„ = 1.75 nm. 



create highly-diffracting conditions and thus appear dark in 
the bright-field TEM images. Nearby grains misoriented by 
only a few degrees are not in the highly diffracting conditions 
and appear light in contrasts. Of particular note is the grain 
boundary grooving observed between these grains. The grain 
boundary morphology or grooving is used to define grain sizes 
in the other parts of the image where adjacent grains are not 
distinguishable by diffraction contrast. From this contrast and 
microstructure, we infer grain sizes in this image to be on 
the order of 15-20 nm. The Z-contrast STEM image of the 
Fig. [TJc) shows the lattice fringes of the NbN film and the 
ability of the film to conform to surface irregularities on the 
substrate. The high-resolution TEM image of Fig. [TJd) was 
taken from a sister sample to the sample shown in Fig. [TJa)- 
(c). In this image, it appears that the top layers of the substrate 
have been modified by the pre-nitridation step used in the film 



preparation. Laterally, the film interface is abrupt and well- 
defined. Grain boundaries can again be identified by changes 
in the atomic structure and the grain boundary grooving de- 
scribed above. Grain sizes in this image are in the 5-20 nm 
range and are typical of the films grown in this series. The 
grain boundary structures shown here are important for pro- 
viding the collective pinning of the vortex lattice in the super- 
conducting state and determine the value of the critical cur- 
rent, as will be revealed in the next section in the analysis of 
the IV characteristics. 

The grooving along grain boundaries can provide sufficient 
contrast in high magnification, plane-view SEM images to en- 
able a much broader analysis of grain sizes in these films. An 
in-depth analysis of the large 1 fim 2 SEM image in Fig. |TJa) 
is performed after converting the gray-scale image of differ- 
ent grain orientations to a binary black and white image in 
Fig. |2ja). This allows a better delineation of grains (black 
domains) for further identification. Because the identification 
of domains depends on the specific threshold used for creat- 
ing a black and white image one should keep in mind that 
our quantitative analysis provides only an approximate count 
of domains that in principle should be verified by a series of 
TEM cross-sectional images. For a quantitative analysis of 
domains in the large 1 fj.m 2 view area, we used the "Analyze 
Particles" method of the image tool Imagefi^ 2 ^ to generate a 
histogram of grain size area A of the black domains with bins 
of size 10 nm 2 . The corresponding linear size, L - yA, his- 
togram is shown in Fig. |2jb). L varies over several orders of 
magnitude between 1 .75 nm and 425 nm with a mean value of 
Lave = 9.2 nm. The grains that were mapped onto black do- 
mains in Fig. [2ja) account for roughly 47% of the completely 
covered film. 

The key result of the quantitative SEM analysis is that the 
distribution of grains is dominated by small grains, which are 
reasonably well described by an inverse power-law 1/L 3 , as 
shown by the red solid line in Fig. [2jb). Notably the major- 
ity of grains extracted from the large view SEM image is less 
than 20 nm, which is consistent with the TEM cross-sections 
shown in Fig. |TJb) and (d). The importance of the peculiar 
inverse power-law will be revisited in the analysis of the de- 
pinning current of the flux lattice state in Sec. III.B. 



C. Resistivity 

Standard film characterization was performed by transport 
measurements in the thin film. Figure [3] shows the resistivity 
data close to the superconducting transition T c = 13.7 ± 0.2 
K at zero-applied magnetic field B = 0T. The resistivity is in 
agreement with reports for thicker films in Ref . [28] The tran- 
sition is suppressed with the field applied perpendicular to the 
film. At B = 1 T and 5 T we observe T C (B) = 13.2 K and 
1 1 .2 K, respectively. From this field dependence we estimate 
for the slope of the upper critical field at the phase transition, 
dB C 2/dT - -2.0 T/K, and derive a zero-temperature coher- 
ence length £ = V _( 5o / (2nT c dB C 2 /dT) = 3.5 nm, where Oq 
is the flux quantum. In addition, in superconductors where 
the mean-free-path ( of electrons is shorter than the zero- 
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FIG. 3. (color online). Resistivity p vs. temperature T in magnetic 

fields B = 0, 1,5 T. Inset: p(T) between room temperature and T c is p IG 4 (color online), (a): IV curves at several representative 

characteristic of a bad metal with a RRR value less than one, RRR = magnetic fields. Three different regions in the IV characteristics can 

p(300 K)/p(T c ) = 0.83. be clearly seen in (b). 



temperature coherence length i.e., superconductors in the 
dirty limit, the diffusion constant of electrons can be obtained 
directly from D = -4k B /(ne dB c ildT) = 0.55 cm 2 /s, where 
ks is the Boltzmann constant and e is the electron charge. 

We use electronic structure calculations of the electronic 
dispersion along high-symmetry directions to estimate the 
Fermi velocity in bulk NbN to be of the order vj « 100 - 150 
km/sJ 29 " 32 ! Then from the diffusion coefficient D = Vf-C/3 = 
0.55 cm 2 /s, we obtain for the mean-free-path I w 1.0 - 
1 .5 nm < which is roughly three to ten times larger than 
reports for disordered thick films (d > 50 nm) grown by 
magnetron sputtering. 33 We can perform a consistency check 
to see whether these values compare reasonably well with a 
rough estimation from superconducting parameters by em- 
ploying Vf ~ 7r£oA/7z. If we assume V/ ~ 150 km/s, the super- 
conducting gap A « 3 meVj 34 " 36 ! and the relation between the 
clean and dirty limit coherence length in superconductors^ 
£o — l-4£ 2 /f, then for £o = 3£ the mean-free-path is i m 
0.46 £ = 1.6 nm. In conclusion, all these estimates are in 
agreement with each other and our earlier analysis of the film 
morphology of ENABLE grown NbN films, which points to- 
ward a superconductor in the dirty limit. 

The inset in Fig.[3]shows the temperature dependence of the 
resistivity p(T) up to 300 K. We extract a resistivity ratio of 
RRR=0.83 between room temperature and slightly above T c , 
which is indicative of charge transport in a bad metal, where 
defect and grain boundary scattering are important. Such a 
scenario is consistent with a short mean-free-path I < £o and 
superconductivity in the dirty regime. 



III. TRANSPORT MEASUREMENTS IN THE 
FLUX-FLOW STATE 

In type-II superconductors (A » £) vortices enter the super- 
conductor above the lower critical field B c \. Under an applied 



DC transport current, vortices are driven by the Lorentz force 
perpendicular to the current, which is balanced by the pinning 
force in inhomogeneous superconductors as long as the trans- 
port current is smaller than the depinning current. In this case, 
vortices do not order into a vortex lattice state due to pinning 
at inhomogeneities. Finally, when the transport current ex- 
ceeds the depinning current, vortices move, resulting in the 
flux-flow state. In the flux-flow state, the inhomogeneities in 
the superconductor are quickly averaged out by vortex motion 
and the lattice order is recovered! 38 ! 39 ! 

In our experiments, sample environments at fixed temper- 
ature down to 2.3 K were controlled in a Janis flow cryostat 
in low DC magnetic fields, from to 0. 1 T, generated by an 
electromagnet at room temperature and in a single shot fridge 
in an Oxford superconducting magnet with a Variable Tem- 
perature Insert for fields between 0.1 T and 15 T. The critical 
current measurements were carried out by means of two dif- 
ferent pulsed techniques to avoid damage, self-heating and/or 
thermal runaway of the samples. We varied pulse duration and 
duty cycle as a mean of assessing and minimizing the self- 
heating. In the first method the commercial set of Keithley 
instruments, Nanovoltmeter Model 2 128 A and Model 6221 
AC/DC Current Source, was used in synchrony to create the 
characteristic IV curves of the superconducting NbN wires at 
fixed magnetic fields and temperatures. In the second method 
two waveform generators were used to create a periodic signal 
applied to the sample through a shunt resistor used to record 
the current. Voltage was measured from the corresponding 
sample leads, both signals were linearly amplified and stored 
in a scope. Both pulsed techniques show results in excellent 
agreement. The dimensions of the NbN film are I (length) x 
w (width) x d (thickness) = 100 pirn x 5 fim x 11.2 nm. All 
IV transport measurements were performed at T — 2.3 K. For 
our analysis we use £ « 3.5 nm. In addition, from the London 
penetration depth, A « 410 nm, reported for films of similar 
thicknessj 3 ^ we deduce the Pearl length A = 2A 2 /d * 30 yt/m 
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FIG. 5. (color online). Dependence of the critical velocity on 
magnetic field in the LO instability region. In the low-field region 
v lo ~ 1 / Vfl, whereas in the high-field region v LO is independent of 
B. The solid lines are guide to the eyes. 



relevant for screening of magnetic flux in thin films. Typi- 
cal IV curves at several magnetic fields are depicted in Fig. [4] 
Three different regions can be clearly seen. Below the depin- 
ning current Id vortices are pinned and the superconductor is 
in the zero-voltage state. Above the depinning current, the 
vortex lattice moves in the flux-flow state causing dissipation. 
At a critical current 7 L o (voltage V L o) an instability occurs and 
the superconductor switches to the normal state. The depen- 
dence of Id and 7lo on the external magnetic field is nonlinear, 
and both decrease with field as - ln(B) as will be discussed be- 
low. As B — > 0, Id and I L o become the same and they are close 
to the departing current. 



A. Larkin-Ovchinnikov instability 



The jump at I L o is due to the instability of the collective mo- 
tion of the vortex lattice as predicted by Larkin and Ovchin- 
nikov (LO) several decades ago. 40 The instability is related to 
the quasiparticle relaxation, thus one can extract the inelastic 
quasiparticle relaxation time in the magnetic field from the in- 
stability. The argument for this effect is as follows: as vortices 
move, an electric field is induced in the normal core of vortices 
resulting in dissipation. Additionally, the electric field shifts 
the distribution of quasiparticles and pushes them outside the 
normal core. As a consequence, the size of the vortex core 
shrinks 



1 + v 2 /v 2 
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where f (0) is the coherence length at the velocity v = and 
v L o is the critical velocity of the instability. Because of the 
reduction of the size of the vortex core, the Bardeen-Stephen 
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For a given Lorentz force, the increase of vortex velocity di- 
minishes T](v), and hence increases v even further. This pos- 
itive feedback speeds up the vortex motion and renders the 
flux-flow state unstable at a critical velocity 
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l LO 



= 1.31- Vl -T/T c . 



(3) 



where D = v/(/3 is the quasiparticle (electron) diffusion con- 
stant with Fermi velocity v/. Thus one can determine t by 
measuring the instability velocity of the flux-flow state. The 
LO instability of the flux-flow state has been observed both i n 
conventional 4 ^" 4 ^ and high-r e cuprate superconductors EES 

Typical IV curves at several magnetic fields are depicted in 
Fig. [4] Above the depinning current Id, vortices move giving 
rise to the flux-flow state. The resistance in the flux-flow state 
increases with current because of the shrinkage of the vortex 
core as a result of the nonequilibrium LO effect . At a critical 
current (voltage), the system switches to the normal state. The 
critical velocity vlo is given by vlo = Vlo/(jj.qBI) with Vlo 
the voltage at the end of the flux-flow branch, B the applied 
magnetic field, and fiQ the vacuum permeability. At higher 
field, the transition to the normal state becomes smooth due 
to the increase of vortex viscosity. We extract Vlo from the 
IV characteristics using the criteria that at V L o the derivative 
dV/dl jumps. 

The dependence of vlo on B is shown in Fig. [5] It de- 
creases as 1 / VB for a weak B and then saturates at a constant 
value. By assuming a uniform distribution of quasiparticles in 
superconductors, LO predicted that the critical velocity vlo is 
independent on the applied magnetic field, see Eq. ([3]). The 
uniform distribution of quasiparticles is realized at high mag- 
netic fields where the inter-vortex distance a is small, such 
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FIG. 6. (color online). Extracted inelastic quasiparticle relaxation 
time r is independent of B in the high field region. 
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second term accounts for the pinning effect, and the last term 
is the damping due to the suppression of superconductivity 
around the vortex core. Here R n w 2.8 kQ is the normal-state 
resistance at T c and a, B, c are fit parameters that depend on B. 
In the limit V <K VLo> we should recover the linear IV curve 
for I— Id = V/Rf where Rf = R„B/B C 2 is the Bardeen-Stephen 
flux-flow resistance. We thus expect c — 1 and a(B) = B/B C 2- 
Our experimental data can be fitted very well by Eq. Q as 
presented in Fig. [7] The fit parameter a{B) is shown in Fig. [8] 
where a(B) * B^/B as expected for high fields not too close 
to B C 2 and from the limit of V/Vlo ^ 1- The other fit pa- 
rameter c x 1 . From the fitted curves, an unstable branch of 
the IV curve near the LO instability becomes visible and the 
system develops hysteresis around the instability region. The 
hysteretic IV curve, due to the LO instability, has also been 
observed experimentally in Ref.l42l 



FIG. 7. (color online). LO instability fits (lines) according to Eq. Q 
of experimental (symbols) VI curves. 



B. Vortex pinning at grain boundaries 



that vloT » a with a * -\/(£> /B. In the low magnetic field re- 
gion, Eq. ^ becomes inapplicable, because the quasiparticle 
distribution is nonuniform and is confined inside the unit cell 
of the vorltex lattice at v* Q given by Eq. d3l, i.e. v* r <K a. 
When the velocity of vortices increases such that quasiparti- 
cles are no longer confined in the unit cell of vortex lattice, 
i.e. the condition vr » a is fulfilled, the flux-flow instabil- 
ity is triggered 49 . Therefore vlo ~ ci/t in the low magnetic 
field region, as shown in Fig. [5] A magnetic-field indepen- 
dent critical velocity Vlo m the high-field region indicates that 
the heating effect due to vortex motion is weak and can be ne- 
glected. In the opposite case of large self-heating, it was found 
theoretically 50 and experimentally 5 -^ that vlo decreases as 
the magnetic field increases, vlo ~ 1/ Vs, which is clearly 
not the case here. 

In the next step, we use Eq. Q to find the inelastic relax- 
ation time for quasiparticles at high fields, as shown in Fig. [6] 
The plateau of t ~ 2 ns above 8 T is expected for flux-flow 
dominated by the LO instability at high fields, while the ori- 
gin of the rise in t between 6 T and 8 T is not understood at 
this time. A similar plateau in the relaxation time of NbN 
was reported in Ref. 46 Obviously, for technical applica- 
tions a shorter relaxation time is preferred in order to achieve 
a faster response of the NbN SNSPD after the formation of 
a hot spot. For that reason the authors of Ref. |46] fabricated 
NbN/ferromagnetic hybrids, where t is two orders of magni- 
tude smaller than in conventional NbN film due to the addi- 
tional scattering channel of quasiparticles by magnetic impu- 
rities in the ferromagnetic layer, leading to a faster relaxation 
of excited nonequilibrium quasiparticles. 

The TV curves including the LO instability can be described 
by the following phenomenological equation 44 



/(V) = 



V 

Rn 



a(B) 



< B(B)(V/V LO y 

1 + v z /v 2 T 1 + v 2 /v'. 



+ 1 



LO 



(4) 



where the first term in the bracket accounts for the reduced 
Bardeen-Stephen viscosity in the nonequilibrium region, the 



We proceed with the investigation of the depinning transi- 
tion of vortices. We use the practical criterion that a depin- 
ning transition occurs when the measured voltage V is larger 
than an arbitrary threshold of 1 mV. From that we can obtain 
the dependence of the depinning current Id on magnetic field 
from the IV curves. The results are compiled in Fig. [9] The 
depinning current Id depends weakly on B when B* < 10 mT 

and decreases as Id In B above. It is worth to note that this 

logarithmic dependence cannot be explained by the collective 
pinning theory,^ which predicts Id ~ B 2 for thin films. 

It is known from thin-film growth parameters and con- 
firmed by HTEM and SEM images that our NbN films ex- 
hibit island-like growth resulting in granular morphology. 
Since superconductivity is suppressed at grain boundaries, 
they provide a pinning potential for vortices and thus may af- 
fect SNSPD detector performance. The pinning due to grain 
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FIG. 8. (color online). Dependence of a(B) on B obtained from fits 
of Eq. Q to experimental data. The solid line is a guide-to-the eye 
for the Bardeen-Stephen flux-flow behavior for B c \ «B< B C 2- 
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FIG. 9. (color online). Experimental depinning current lj as a func- 
tion of B (symbols) at low DC magnetic fields (a) and high magnetic 
fields (b). The high field data are measured in a single-shot fridge 
in an Oxford superconducting magnet. Lines are theoretical value 
obtained by using Eq. l|6l and the distribution function in Fig. [10] 



boundaries has also been observed experimentally in high- 
temperature superconductor YBa2Cu307_^ filmsj^ where the 
dependence of Id on B is similar to that in Fig. [9] The 
depinning current density at B = T of our thin film is 
jd - Id/(wd) - 7.7 MA/cm 2 , which is close to the depair- 
ing current density j c = cO /(12^ 2 V3/l 2 f) « 16 MA/cm 2 . 
Thus the pinning of vortices by grain boundaries can achieve 
a depinning current that is close to the depairing current. 

The depinning current depends crucially on the size dis- 
tribution of superconducting grains (domains), which can be 
extracted from the dependence of Id on B. We use a model for 
collective breakaway of pinned vortices developed by Fedotov 
et For simplicity we further assume square domains in 
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FIG. 10. (color online). Size distribution of grains obtained from Eq. 
|7](. Inset: Semi-logarithmic plot of the same distribution. 



our thin-film NbN. A pinning theory for more general shapes 
of domains was presented in Ref.[55]and the results are quali- 
tatively similar to those with square domains. One starts with 
a given probability density distribution P(L) of grains with 
linear dimension L, where the probability density for a vortex 
inside the square domain of size L is W(L) = NL 2 P(L) with 
normalization constant N. The probability for finding a sin- 
gle vortex j>inned in an infinitely large domain is normalized 
to unity, f Q W(L)dL= 1. 

Not all vortices can sit at the energy minimum of the 
pinning potential due to the competition between the pin- 
ing energy and elastic energy, and the resulting vortex con- 
figuration is a compromise between these two energies. 
Assuming a square pinning potential with strength e p = 
O 2 ) r 2 /(87T/I^o) 2 an d using the expression for the elastic energy 
e e = <f>oBd 2 /(8777I) 2 , we obtain the maximal displacement for 
pinned vortex, 6(B) - ^nr^H^JB, i.e., if the displacement 
of the vortex from the grain boundary is less than 6, then the 
vortex remains pinned. Here r c ~ £0 characterizes the strength 
of pinning potential. Under these conditions the vortex core 
gains condensation energy over the elastic deformation energy 
by staying at the grain boundary. It follows that the probabil- 
ity of a vortex to lie less than a distance +6 away from the 
grain boundary is equal to the ratio of the area of each strip to 
the total area 



L 2 



if L < 26, 
if L > 26. 



(5) 



For this case it was shown that the depinning current, normal- 
ized to its value at B = 0, is simply the ratio of pinned vortices 
n p to the total number of vortices n, nt , 



Id(B) 
Id(0) 



n,ot Jo 

T 

J 25 



W(L)P(L; 6)dL 



1 - AM P(L)(L - 26fdL. (6) 
lis 

Here the magnetic field enters only through the vortex dis- 
placement 6(B). 

Since a weak magnetic field corresponds to soft elastic 
shear stiffness of the vortex lattice, i.e., 6 ~ yf^jAB, low- 
field measurements probe primarily large domains. From the 
measured Id as a function B, we can obtain the distribution 
function from Eq. (|6]l 



P(y) 



1 



d 3 i d (y) 



2NI d (0) dy 3 



(7) 



with y = 26. To get rid off the small oscillations in Fig. [9] we 
first smooth the experimental data and then use Eq. (j7]). The 
resulting distribution function is shown in Fig. [10] which is 
in reasonable agreement with the experimental one shown in 
Fig. [2] The probability distribution for small grains is much 
larger than that for larger grains. The distribution function for 
small grains follows 1/L 3 , while for large grains, it follows 
an exponential distribution. Since the smallest length scale 
for vortex is of order of we cannot resolve the distribution 
for grain sizes smaller than £ from the measurements in the 
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flux-flow state. We confirm the agreement between the mea- 
sured and calculated critical current by inserting the extracted 
distribution of Fig. 10 into Eq. d3l, as shown in Figs. [9] 



IV. SUMMARY 

In summary, we characterized the quality and uniformity 
of ENABLE-grown thin-film NbN superconductors for poten- 
tial SNSPD applications. From transport measurements we 
derived superconducting material parameters T c - 13.7 K, 
£, — 3.5 nm, and depinning current density = 7.7 MA/cm 2 . 
In addition, we determined that our thin films of thickness 
d = 11.2 nm are in the dirty limit with the mean-free-path 
much shorter than the coherence length of the hypothetically 
clean superconductor, I w 0.15fo- This length scale was fur- 
ther corroborated by the distribution of grain sizes extracted 
from the analysis of an SEM image and the field-dependence 
of depinning currents. The presented vortex theory success- 
fully explained the collective vortex lattice motion in the flux- 
flow state with the Larkin-Ovchinnikov instability in the IV 
characteristics at high bias currents, as well as the depinning 
current 1 j at low bias currents. The detailed analysis of the 
Larkin-Ovchinnikov instability revealed a relatively long in- 
elastic quasiparticle relaxation lifetime of order ~ 2 ns, which 
might provide the bottleneck for the hot spot relaxation in 
NbN-based SNSPD devices. Finally, from the nondestructive 



measurement of the depinning current with magnetic fields, 
we extracted the characteristic domain size distribution of 
grains, which resulted in comparable values to the indepen- 
dent, yet destructive, analysis using HTEM techniques. While 
the prevalence of grain boundaries in thin-film NbN supercon- 
ductors crucially affects transport properties like critical cur- 
rents, their potential for vortex pinning at low bias currents 
is negligible for SNSPD applications, which are typically bi- 
ased close to the critical current. This study has shown the 
potential use of field-dependent measurements of the depin- 
ning current in micron-sized wires for determining the grain 
size distribution in thin-film superconductors. The advantage 
of a nondestructive characterization method of the uniformity 
of thin superconducting films may prove beneficial for the pre- 
screening of films for further nano-patterning. 
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